Astronomy &: Astrophysics manuscript no. shuliu'17909 
March 28, 2012 



© ESO 2012 



(N 

O 

(N 



(N 

6 

in: 



> 

(N 

oo 

rn 

O 



% 



Abundances of neutron-capture elements in G 24-25 * 

A halo-population CH subgiant 

S. Liu^'2, P. E. Nissen^, W. J. Schuster^ G. Zhao^ Y. Q. Chen^ and Y. C. Liang^ 

^ Key Laboratory of Optical Astronomy, National Astronomical Observatories, CAS, 20A Datun Road, Cliaoyang 

District, 100012, Beijing, China 

e-mail: gzhaoOnao . cas . en 
^ Graduate University of the Chinese Academy of Sciences, 19A Yuquan Road, Shijingshan District, 100049, Beijing, 

China 
^ Department of Physics and Astronomy, University of Aarhus, DK-8000 Aarhus C, Denmark 

e-mail: penOphys .au.dk 
* Observatorio Astronomico National, Universidad National Autonoma de Mexico, Apartado Postal 877, C.P. 22800 

Ensenada, B.C., Mexico 

e-mail: schuster@astrosen.unam.mx 

Received 18 August 2011 / Accepted 5 March 2012 

ABSTRACT 

Aims. The differences between the neutron-capture element abundances of halo stars are important to our understanding 
of the nucleosynthesis of elements heavier than the iron group. We present a detailed abundance analysis of carbon 
and twelve neutron-capture elements from Sr up to Pb for a peculiar halo star G 24-25 with [Fe/H] = —1.4 in order to 
probe its origin. 

Methods. The equivalent widths of unblended lines are measured from high resolution NOT/FIES spectra and used to 
derive abundances based on Kurucz model atmospheres. In the case of CH, Pr, Eu, Gd, and Pb lines, the abundances 
are derived by fitting synthetic profiles to the observed spectra. Abundance analyses are performed both relative to the 
Sun and to a normal halo star G 16-20 that has similar stellar parameters as G 24-25. 

Results. We find that G 24-25 is a halo subgiant star with an unseen component. It has large overabundances of carbon 
and heavy s-process elements and mild overabundances of Eu and light s-process elements. This abundance distribution 
is consistent with that of a typical CH giant. The abundance pattern can be explained by mass transfer from a former 
asymptotic giant branch component, which is now a white dwarf. 
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1. Introduction 

The heaviest elements in halo stars are mainly produced by 
neutron captures divided into the slow (s-) and the rapid 
(r-) processes, which have different reaction timescales. 
Among the s-process synthesized elements, lead (Pb) is 
particularly interesting because it is the most abundant cl- 
ement in the third peak of the nuclei distribution in stars. 
However, Pb lines are very weak in low metallicity stars, 
and thus Pb abundances have not been measured for many 
stars compared with other s-process elements, such as Ba. 
Pb is often detected for carbon-enhanced stars with en- 
hanced Ba abundances. For example, Bisterzo ct al. (2006) 
compiled a catalog of 23 Pb-rich stars found in the litera- 
ture. Most of these stars have [Fe/H] < —2 and are carbon- 
enhanced metal-poor stars with s-process enhanced abun- 
dances (CEMP-s) or with r- and s-process enhanced abun- 
dances (CEMP-r/s). Eleven stars of the CEMP-s/CEMP- 
r/s compiled in Bistcrzo's work are proven members of 
binary systems (McClurc & Woodsworth 1990; Aoki et 



Based on observations made with the Nordic Optical 
Telescope on La Palma. Table 2 is available in electronic form 
at http : //www . aanda . org 



al. 2001, 2002b; LucatcUo et al. 2003; Cohen ct al. 2003; 
Sivarani et al. 2004) . 

In the past ten years, some studies also reported Pb 
measurements for Ba stars and CH giants and subgiants 
(Van Eck et al. 2001, 2003; Johnson & Bohe 2004; Allen 
& Barbuy 2006; Pereira & Drake 2009, 2011; Goswami & 
Aoki 2010). These classes of stars exhibit overabundances 
of s-proccss elements; in addition, the CH stars have high 
C abundances. Furthermore, these stars fall into different 
metallicity ranges and belong to different populations in 
the Galaxy. Most of the Ba stars are found in the disk 
with [Fe/H] from -0.5 to +0.4; fewer than 6% of the Ba 
stars belong to the halo (Mcnnessier et al. 1997; Allen & 
Barbuy 2006). The classical CH giants are clearly mem- 
bers of the halo population with —1.5 <[Fe/H]< —1.0 
(Vanture 1992a). The CH subgiants, i.e stars with similar 
abundance distributions as the CH giants, but situated on 
either the main sequence or the subgiant branch, have typ- 
ical [Fe/H] values in the range from —0.5 to —0.2 and disk 
population kinematics (Luck & Bond 1982, 1991; Sneden 
& Parthasarathy 1983; Smith ct al. 1993). 

Systematic spectroscopic studies indicate that all Ba 
and CH stars belong to binary systems (McClure & 
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Table 1. Atmospheric parameters and radial velocities of G 24- 
25 and G 16-20. 

Star Tcft logs [Fe/H] ^ ''rad 
(K) (kms~^) (kms~^) 

G 24-25 5828 3.86 -1.40 1.20 -312.9 

G 16-20 5625 3.64 -1.42 1.50 170.8 



Woodsworth 1990). At low metallicity, [Fc/H]< -1.8, 
Lucatello et al. (2005) found that 68% of 19 CEMP-s stars 
display evidence of radial velocity variations. This high 
frequency suggests that all CEMP-s stars are metal-poor 
analogs of the classical CH giants. The origin of the en- 
hanced carbon and s-process elements (including Pb) in 
these stars could be explained by mass transfer from a 
former asymptotic giant branch (AGB) companion (now 
a white dwarf) in a binary system. 

In a study of two distinct halo populations in the solar 
neighborhood, Nissen & Schuster (2010, 2011) measured 
the abundance ratios and kinematics of 94 dwarf and sub- 
giant stars and found a strong enhancement of [Ba/Fe] for 
a halo star, G 24-25, with [Fc/H] = —1.4. Further investiga- 
tion of the spectrum identified a Pb line at 4057.8 A, which 
is not seen for the rest of the stars. For example, the Pb fea- 
ture at 4057.8 A cannot be detected in G 16-20, a halo star 
with a similar metallicity as G 24-25. In addition, G 24-25 
has a strong CH band near 4300 A, indicating that it is a 
CH subgiant. To our knowledge, this is the first known CH 
subgiant with a metallicity that is typical of the halo pop- 
ulation. In this paper, we report our abundance measure- 
ments of Pb and other neutron-capture elements for this 
peculiar halo star. For comparison, G 16-20 is analyzed as 
a reference star. The abundance pattern is compared with 
other chemically peculiar stars in the literature, and simple 
AGB wind- accretion models are adopted to reproduce the 
measured abundances. 



2. Abundance analysis 

2.1. Observational data and atmospheric parameters 

As described in Nissen & Schuster (2010), the spectra of 
G 24-25 and G 16-20 were obtained at the 2.56m Nordic 
Optical Telescope (NOT) using the Fibre fed Echelle 
Spectrograph (FIES) with a resolving power of R ~ 40,000. 
The spectral range goes from 4000 A to 7000 A, and the 
signal-to-noise ratio (S/N) equals about 170 at 5500 A but 
is only between 50 and 80 in the blue spectral region. 

The stellar atmospheric parameters for G 24-25 and 
G 16-20 are taken from Nissen & Schuster (2010) and given 
in Table 1, where the heliocentric radial velocities are also 
given. According to Nissen & Schuster (2010), the total 
space velocities with respect to the local standard of rest 
are 315 km s'^ and 263 km s'^ for G 24-25 and G 16-20, re- 
spectively, which clearly classify both stars as members of 
the halo population. 

2.2. Line selection and atomic data 

The atomic lines were selected from previous abundance de- 
terminations of the heavy elements of stars, namely Sneden 
& Parthasarathy (1983), Gratton & Sneden (1994), Sneden 
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Fig. 1. The CH feature near 4310 A. The observed spectrum of 
G 24-25 is shown as bold dots. The solid line shows the synthetic 
spectrum corresponding to [C/Fe] = 1.05. The dot-dashed and 
dashed lines, corresponding to A[C/Fe] = ±0.07, are shown to 
illustrate the sensitivity of line strengths to variations in the 
carbon abundance. 



et al. (1996), Reddy et al. (1997), Aoki et al. (2001, 2002b), 
Cowan et al. (2002), Sivarani et al. (2004), and Johnson 
& Bolte (2004). The log g/ values are primarily adopted 
from high-precision laboratory measurements; references 
are given in Table 2. 



2.3. Abundance calculations 

For most elements, the abundances are derived from equiv- 
alent widths (EWs) of unblended lines measured by fitting 
a Gaussian to the line profile. For blended lines and/or lines 
with significant hyperfine structure (HFS) and/or isotope 
splitting, the abundances are derived by spectrum synthe- 
sis. The model atmospheres of Kurucz (1993) are used, and 
local thermodynamic equilibrium (LTE) is assumed. 

The carbon abundance is derived from the EWs of 
three unblended Cl atomic lines providing [C/Fe] = 1.03 
and [C/Fe] = -0.39 for G 24-25 and G 16-20, respectively. 
Compared to G 16-20, G 24-25 has clear carbon atomic lines 
and strong CH and C2 lines. From spectrum synthesis fit- 
ting of the CH A-X band of G 24-25 around 4310 A (Fig. 1), 
we obtained [C/Fe] = 1.05 when using the CH molecular 
line data of Barklcm et al. (2005). The abundances derived 
from the C i atomic lines and the CH A-X band are con- 
sistent, whereas Sivarani et al. (2004) reported that carbon 
abundances derived from C i lines are on average about 0.1- 
0.3 dex higher than those obtained from the CH lines. 

Since the CN 3800 A band is not covered in our spec- 
tra, we tried to derive the nitrogen abundance by spectrum 
synthesis fitting of the CN 4215 A band adopting molecu- 
lar line data from Aoki et al. (2002a). The CN lines around 
4215 A are quite weak (EW < 5 mA), and due to the rather 
low S/N ratio in this spectral region, we can only get an 
upper hmit [N/Fc] < 0.16 dex for G 24-25. In the case of 
G 16-20 we get [N/Fc] < 0.0 dex. 

The abundances of the light s-process elements Sr^, Y, 
and Zr, as well as the heavy neutron-capture elements Ba, 



^ We do not include the Sril line at 4161.8 A, because the 
abundance inferred from that line is 0.3 dex higher than the 
average for the other lines (Sneden et al. 1996; Aoki et al. 2001). 
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La, Ce, Nd, and Sm arc derived from EW measurements, 
whereas spectrum synthesis is used for Pr, Eu, Gd, and Pb. 
The HFS of the Ban lines at 5853.7, 6141.7, and 6496.9 A 
is not taken into account since it is insignificant according 
to Sneden et al. (1996) and Mashonkina et al. (2010). Our 
Ba abundances are in good agreement with those in Nissen 
& Schuster (2011), who also found that the effects of HFS 
on the Ba lines are negligible. The La ii and Sm ii lines are 
quite weak, and we do not include any HFS. 

The collisional broadening of lines induced by neutral 
hydrogen is also considered. The width cross-sections of 
the Ci, Sri, Srii, and Ban lines are taken from Anstee 
& O'Mara (1995), Barklem & O'Mara (1997, 2000), and 
Barklem et al. (2000). For the remaining lines, we follow 
Cohen et al. (2003) and adopt the Unsold (1955) approxi- 
mation to the van der Waals interaction constant enhanced 
by a factor of two. For the four strong lines, Sm 4077.7 A 
and the Ban lines at 5853.7, 6141.7, and 6496.9 A, the 
effect on the derived abundances of the uncertainties in 
the damping constants is approximately 0.08 dex. However, 
most element abundances in our work are based on weak 
lines {EW < 60 mA), for which the influence on the results 
caused by the uncertainties in the collisional cross-sections 
is negligible. Hence, we do not discuss possible errors in the 
damping constants in Sect. 2.4. 

The spectrum synthesis fits to the Pr, Eu, and Gd 
lines are shown in Fig. 2. The atomic and hypcrfine struc- 
ture data for the Pri 5322.8 A and Eun 4129.7 A lines are 
adopted from Mashonkina et al. (2009) and Mashonkina & 
Gehren (2000), respectively. The log g/ value for the Gdn 
line at 4191.1 A is taken from Den Hartog et al. (2006). 

In the observed spectrum of G 24-25, we could detect 
one lead line, the Pb i A4057.8; which is too weak to be seen 
in the spectrum of G 16-20 (Fig. 3). The log 5/ value is taken 
from Biemont et al. (2000) and the hyperfine data are taken 
from Aoki et al. (2001). Compared to G 16-20, a more nor- 
mal neutron-capture star, it is clear that G 24-25 exhibits a 
strong Pb line. The spectral region from 4057.6 A to 4058 A 
does not show a smooth line profile for G 24-25, but con- 
sidering that the estimated S/N around the Pb line is only 
about 60, we propose that this is probably due to noise. 
Thus, a x^ fitting of the ten observed data points from 
4057.69 A to 4057.92 A was applied to determine the lead 
abundance. The x^ ^^^ computed as described by Nissen 
et al. (1999), i.e. 



X 



E(Oi-so^ 



where Oi is the observed relative flux. Si the synthetic 
flux, and a = {S/N)''^ = 1/60. The value of [Pb/Fe] 
was varied in steps of 0.02 dex to flnd the lowest x^- 
The result is a parabolic variation in x^ ^^ shown 
in Fig. 4. The most probable value of [Pb/Fe] corre- 
sponds to the minimum of x^, and Ax^ = 1, 4, and 
9 (the dashed horizontal lines) correspond to the 1-, 
2-, and 3-cr confidence limits of determining exclusively 
[Pb/Fe]. The spectrum synthesis fit to the Pb line for the 
abundance corresponding to the minimum x^ ([Pb/Fe] = 
1.68) is shown as the solid line in Fig. 3. A Gaussian broad- 
ening function was applied to fit the observed spectrum in 
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Fig. 2. Comparison of observed (bold dots) and synthetic pro- 
files of the spectral lines of the heavy elements Pr (top panel), 
Eu (middle), and Gd (bottom). The respective elemental abun- 
dances corresponding to the full drawn lines are listed in Table 5. 
Two alternative synthetic spectra are shown to demonstrate the 
sensitivity of the line strength to the abundance. 



order to take the combined effect of instrumental, macro- 
turbulent, and rotational broadening into account. The 
nearby spectral lines could be well-fitted in this way leading 
to a FWHM of the Gaussian function of 7.8±lkms-i. The 
corresponding uncertainty in [Pb/Fe] is around 0.04 dex. 

2.4. Error analysis 

The uncertainties in the abundances arise from random and 
systematic errors. The random errors associated with the 
errors in the log 5/ values and the EW measurements can be 
estimated from the line-to-line scatter in the abundances of 
a given element, based on many lines. The random error in 
the mean abundance is (Jew/ v^: where N is the number 
of lines studied. For the elements that have only one line, we 
use three times the uncertainty given by Cayrel (1988) by 
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Table 3. Abundance uncertainties for neutron-capture ele- 
ments in G 24-25. 
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Fig. 3. Spectrum-synthesis fitting of the Pb I A4057.81 line for 
G 16-20 and G 24-25. 




[Pb/Fe] 

Fig. 4. The x^ variation in the synthetic fit to ten observed data 
points from 4057.69 A to 4057.92 A applied to the determination 
of the lead abundance in G 24-25. 



taking into consideration the error in the continuum recti- 
fication. For instance, Sri 4607.34 A with S/N = 130 leads 
to an error of 0.6 mA by using the formula (Eq. 7) given by 
Cayrel (1988). We then estimate that the error in the EW 
measurement is 1.8mA for Sri 4607.34 A, vifhich results in 
an error of 0.06 dcx in the Sri abundance calculation. 

A significant contribution to the systematic errors 
comes from the uncertainties in the stellar atmospheric 
parameters adopted in the abundance determination. This 
was estimated by varying the model parameters with ATcs 
= 100 K, Alogg = 0.10 dcx, A[Fe/H] == 0.1 dcx, and A^ 
= 0.30 km s~^. The corresponding abundance uncertainties 
are listed in Tables 3 and 4, and the total uncertainty for 
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Table 4. Abundance uncertainties for neutron-capture ele- 
ments in G 16-20. 
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each element is estimated by the quadratic sum of the at- 
mospheric and random errors. 

An additional systematic error in the derived abun- 
dances relative to the solar abundances is introduced by 
our LTE assumption and the use of one-dimensional model 
atmospheres. These errors tend, however, to cancel for the 
abundances of G 24-25 relative to G 16-20, because the two 
stars have similar atmospheric parameters and also compa- 
rable abundances of Fe and the a-capture elements. 

3. Results and discussion 

3.1. Abundance results 

The abundances for 13 elements are presented in Table 5, 
along with the number of lines on which they are based, 
the adopted solar abundances, log e — log (iVx/A^n) + 12, 
and both [X/Fe] and [X/Ba] for each species. We adopted 
the standard solar abundances of Asplund et al. (2009) as 
a reference for the various elements. As seen, all of the 13 
elements studied are enhanced in G 24-25 relative to both 
the Sun and the reference star G 16-20. 

In Fig. 5, we present the abundance patterns of G 24-25 
and G 16-20. The abundances of Na, Mg, Si, Ca, Ti, Cr, 
Ni, Mn, Cu, and Zn are adopted from Nissen & Schuster 



S. Liu et al.: Abundances of neutron-capture elements in G 24-25 
Table 5. Abundances of neutron-capture elements of G 24-25 and G 16-20. 
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(2010, 2011). The two stars have similar behaviors for both 
the a- and iron-peak elements but show obvious differences 
for C and heavy neutron-capture elements. This can be at- 
tributed to the chemical peculiarity of G 24-25. In particu- 
lar, the second-peak s-process elements (Ba, La, Ce, Pr, Nd, 
Sm, Gd) are more enhanced than the first-peak s-process 
elements (Sr, Y, Zr) in G 24-25, which also has a strong 
carbon enhancement. 

Taking into account the most recent research on metal- 
poor stars, Beers & Christlieb (2005) divided metal-poor 
stars into various subclasses. According to their definitions, 
G 24-25, which is both carbon- and barium- enhanced with 
[C/Fe] > -hl.O and [Ba/Fe] > +1.0 along with [Ba/Eu] > 
-1-0.5, can be classified as CEMP-s. However, most CEMP- 
s stars found in the literature have [Fe/H] < —2, which 
is more metal-deficient than G 24-25 with [Fe/H] = —1.4. 
On the other hand, CH stars have higher metallicitics than 
CEMP stars and occur in binary systems. Given that G 24- 
25 is a single-lined spectroscopic binary (Latham et al. 
2002), we suggest that G 24-25 should be classified as a 
CH subgiant, although it has a lower metallicity than the 
classical CH subgiants belonging to the Galactic disk. In 
this context, the C and s-process clement enhancements 
are consistent with mass transfer from a former AGB bi- 
nary companion (McClure & Woodsworth 1990). 



3.2. Comparing the abundance pattern to model predictions 

To investigate whether the abundance pattern of G 24-25 
can be explained by mass transfer from a former AGB com- 
panion in a binary system, the observed abundances are 
compared with the values predicted by theoretical models. 
In this comparison, the ratios of the abundances of G 24-25 
to those of the reference star G 16-20 are used to ensure that 
the effects of Galactic chemical evolution are minimized. 

Low-mass AGB stars (1.3 ^ M/Mq <, 3) are affected 
by recurrent thermal pulses in the He shell leading to cor- 
related enhancements of freshly synthesized ^^C and s- 
process elements in their photospheres (Straniero et al. 
1995; Gallino et al. 1998). In addition, the more massive 
AGB stars undergo hot bottom burning where ^^C is con- 
verted into ^^N. The minimum mass for this process is 
somewhat uncertain, but at the metallicity of G 24-25, cor- 
responding to Z ~ 0.001, Karakas & Lattanzio (2007) es- 
timate it to be around 3.5 Mq (see their Table 1). Below 
this mass limit, an AGB star can have a C-rich atmosphere 
with no enhancement of nitrogen as found in the case of 
G 24-25 ([N/Fe] < 0.16). Hence, the mass of the AGB star 
that is responsible for the abundance anomalies of G 24-25 
is unlikely to be higher than about 3 Mq] a simple consid- 
eration of the standard initial mass function implies that 
this mass is more likely to be around 1.5 Mq. 

Bisterzo et al. (2010) present theoretical predictions of 
an updated low-mass AGB stellar nucleosynthesis model 
at different metallicitics. For a model with a given initial 
AGB mass and metallicity, they provide the abundances of 
all elements from carbon to bismuth in the envelope. They 
apply the analysis to CEMP-s stars in Bisterzo et al. (2011) 
and get a good fit to the observed abundances by varying 
the "'^'^C-pocket profile and dilution factor. Unfortunately, 
the calculated model data published in Bisterzo et al. (2010, 
2011) are limited (only two ^'^C-pocket cases are available 
for the mass and metallicity we need, M^f^ = 1.5Mo and 
[Fe/H] = —1.3), and the model prediction does not closely 
match the abundances of G 24-25. 

Cristallo et al. (2009) also present a homogeneous set of 
calculations of low-mass AGB models at different metallic- 
itics, which they extend to create a database of AGB nu- 
cleosynthesis predictions and yields, called the FRANEC 
Repository of Updated Isotopic Tables & Yield (FRUITY) 
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database (Cristallo ct al. 2011). This database is an online 
interactive interface of predictions for the surface composi- 
tion of AGB stars undergoing the third dredge-up (TDU) 
when choosing different combinations of initial mass and 
metallicity. The first set of 28 AGB models available from 
FRUITY covers the masses 1.5 < M/Mq < 3.0 and metal- 
licities IxlO'^ < Z < 2x10"^. The stellar models of the 
FRUITY database were computed based on the updated 
FRANEC code (ChiefR et al. 1998). Carbon-enhanced opac- 
ity tables were adopted to take into account the effects of 
TDU episodes. The computation ends when the minimum 
envelope mass for TDU occurs. For our purpose, we adopted 
the AGB model computation from FRUITY for an initial 
mass M = 1.5Mq and Z = 0.001 (corresponding to [Fe/H] 
= —1.3) and scaled the whole abundance pattern down to 
our observational data by a factor of 1.6. These AGB model 
predictions for neutron-capture elements from Sr to Pb is 
presented by a solid line in Fig. 6. For most elements, the 
observed abundances are fitted quite well, but the Eu abun- 
dance is significantly overabundant according to the model 
prediction. 

Alternatively, the parametric AGB model for metal- 
poor stars from Zhang et al. (2006) and Cui et al. (2010) 
has been used to predict the abundance pattern of G 24-25, 
although this model is only available for a mass of 3 Mq and 
Z = 0.0001. They adopted the model for metal-poor stars 
presented by Aoki et al. (2001) after updating many of the 
neutron-capture rates according to Bao et al. (2000). This 
approach is not based on detailed stellar evolution models, 
but has been used to successfully explain the abundance 
patterns of some very metal-poor stars. There are three 
parameters in the model: the neutron exposure per ther- 
mal pulse. At, the overlap factor, r, and the component 
coefficient of the s-process, Cg- In the case of multiple sub- 
sequent exposures, the mean neutron exposure is given by 
To = —Ar/hir. 

In their model, the convective He shell and the envelope 
of the giant at some time on the AGB, will be overabun- 
dant in heavy elements by factors of /shell and /cnv,i, respec- 
tively, with respect to Solar System abundances normalized 
to the values of the metallicity, Z. The approximate rela- 
tion between /shell and /env.i is given by Eq.(3) of Zhang et 
al. (2006), where AMdr is the total mass dredged up from 
the He shell into the envelope of the AGB star, and Aff is 
the envelope mass. Given a mass of 3Mq and Z = 0.0001 
for the AGB star, AMdr/Mf is 1/9 according to Karakas 
& Lattanzio (2007). 

For a given s-process element, the overabundance fac- 
tor /env,2 in the companion star envelope can be approxi- 
mately related to the overabundance factor /env,i by Eq.(4) 
of Zhang et al. (2006), where AM2 is the amount of matter 
accreted by the companion star from the AGB progenitor, 
and Af| is the envelope mass of the accreting star. The 
component coefficient, Cs = /env,2//shoii, is then computed 
from Eq.(5) of Zhang et al. (2006). 

The closest possible match between the observed and 
predicted abundances is obtained for the following values 
of the three parameters: overlap factor, r = 0.50, neutron 
exposure. At = 0.26 mb~^ (i.e. a mean neutron exposure 
To = 0.375inb^^), and s-process component coefficient, 
Cs=0. 0012384. The predicted abundance pattern is given 
as the dashed line in Fig. 6. 

To highlight the AGB contribution. Fig. 6 shows 
[X/Fejc 24-25 rninus [X/Fe]Gi6_2o as a measure of the 
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Fig. 6. The observed abundances of the neutron-capture el- 
ements of G 24-25 relative to G 16-20 in comparison with two 
AGB model predictions. 



abundances of the atmospheric material that has been ac- 
creted by G 24-25 from its donor star. We note, however, 
that the values plotted for [Pr/Fe] and [Pb/Fe] are lower 
limits to the difference between G 24-25 and G 16-20, be- 
cause the Pr and Pb lines could not be detected in G 16- 
20, hence we have only upper limits to both [Pr/Fe] and 
[Pb/Fe] for this reference star. From Fig. 6, we can see that 
the observed abundance distribution is consistent with the 
two model predictions, except that a too high Eu abun- 
dance is predicted by the FRUITY model. The consis- 
tency between the observed abundance distribution and the 
model predictions indicates that the overabundances of s- 
process elements in G 24-25 can be assumed to originate 
from mass transfer in a binary system containing material 
produced by an AGB companion. 

3.3. Comparison of G 24-25 with other peculiar stars 

It is also interesting to compare the abundance peculiarities 
of the neutron-capture elements in G 24-25 with those of 
some well-studied CEMP-s or CH stars. We chose from the 
literature stars with as many heavy-element abundances 
as possible and atmospheric parameters close to those of 
G 24-25. Unfortunately, not one chemically peculiar star 
was found to satisfy all these criteria. Instead, two more 
metal-poor stars, HE 0024-2523 (Lucatello et al. 2003) and 
BD-h04°2466 (Pereira & Drake 2009), were considered, 
both of which have 17 elements available for comparison 
and are members of binary systems. HE 0024-2523 is a 
CEMP-s dwarf star with [Fe/H] = -2.7, and BD-h04°2466 
is a classical CH giant with [Fe/H] = —1.9. As seen from 
Fig. 7, G 24-25 shows much the same abundance pattern as 
the two stars, ahhough G 24-25 is closer to BD+04°2466, 
the classic CH star, than to HE 0024-2523. We note that 
HE 0024-2523 was classified as a lead star with an extreme 
carbon enhancement given by [Pb/Fe] = -1-3.3 and [Pb/Ba] 
= +1.9. 



3.4. Comparison of [s/Fe] versus [Fe/H] for chemically 
peculiar binary stars 

Due to the similar formation mechanisms, we have com- 
pared G 24-25 with other classes of chemically peculiar bi- 
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Fig. 7. Comparison of G 24-25 witli HE 0024-2523 (CEMP-s 
dwarf) and BD-F04°2466 (CH giant). 



nary stars, i.e. the CEMP dwarfs, subgiants and giants; the 
CH subgiants and the classical CH giants; the Ba dwarfs 
and giants. The selected binary stars either have derived or- 
bit solutions or exhibit radial velocity variations beyond Scr 
(Lu et al. 1987; McClure & Woodsworth 1990; Preston & 
Sneden 2001; Aoki et al. 2003; Cohen et al. 2003; Lucatello 
et al. 2003; Pourbaix et al. 2004; Tsangarides et al. 2004). 
These various classes are indicated by different symbols 
in Fig. 8, where [hs/Fe] and [Zs/Fe] represent the average 
abundances relative to Fe of the second-peak elements (Ba, 
La, Ce, Nd, Sm) and the first-peak s-process elements (Sr, 
Y, Zr), respectively. Pr and Gd are not included in the 
comparison because most works in the literature have no 
Pr or Gd measurements. We do not include the so called 
"metal-deficient Ba stars" (mdBa). In the past two decades, 
only five mdBa stars have been classified based on chemical 
analyses (Luck & Bond 1991; Junqueira & Percira 2001). 
Some works have presented different views of these five 
"mdBa" stars. BD -^04°2466 was confirmed as a CH giant 
by Pereira & Drake (2009). Jorissen ct al. (2005) considered 
HD 104340 and BD +03°2688 as intrinsic Ba stars, instead 
of classical Ba stars belonging to a binary system, because 
they lie on the thermally pulsing AGB sequence and did 
not show any indication of orbital motion from radial ve- 
locity measurements. Jorissen et al. (2005) also doubt the 
binary nature of HD 55496 and HD 206983 since they had 
no, or inconclusive, radial velocity data. 

For a wide range of [Fe/H] (-3 < [Fc/H] < 0), [Is/Fc] 
has no dependence on [Fe/H] for all classes of chemically 
peculiar binary stars, whereas [hs/Fe] shows a decreasing 
trend with increasing [Fe/H] albeit with some scatter. The 
data of G 24-25 are consisted with this trend indicating 
similar origins for these stars. This decreasing trend with 
[Fe/H] becomes even more obvious in the [hs/ls] versus 
(vs.) [Fc/H] panel. In addition, there is a hint that the 
[hs/ls] {[hs/ls] = [hs/Fc] — [Is/Fe]) ratios of giants are sys- 
tematically higher than those of subgiants/dwarfs at [Fe/H] 
> —0.5, but further checks of the differences in the abun- 
dance analyses between giants and subgiants/dwarfs arc 
needed. 

According to Luck & Bond (1991) and Busso et al. 
(2001), [hs/ls] is an indicator of the neutron-exposure in the 
s-process nucleosynthesis: large [hs/ls] ratios correspond 
to high neutron-exposures. Therefore, the CEMP, CH, and 
Ba stars can be described by the same scenario in differ- 



ent metallicity ranges where the neutron-exposure in these 
stars increases as the metallicity decreases. In agreement 
with this trend. Fig. 9 shows the [C/Fc] ratio as a func- 
tion of metallicity for the same objects as in Fig. 8. All 
chemically peculiar binary stars have a carbon enhance- 
ment that increases with decreasing [Fe/H]. The scatter 
in the [C/Fe] ratios becomes significant for CEMP stars 
with [Fc/H] < —2.0. There is a tendency for the carbon 
enhancement to diverge into two branches as shown by the 
solid (low [C/Fe]) and dashed (high [C/Fe]) lines in Fig. 9. 
Among the high branch, 8 of 11 stars have [Pb/Ba] > 1.0, 
but three stars (LP 625-44, HE 2148-1247, and HD 189711) 
have [Pb/Ba] < 1.0; notably, HD 189711 has large uncer- 
tainties of 0.5 dex in [Pb/Fe] (Van Eck ct al. 2003), and 
HE 2148-1247 is a typical CFMP-r/s star (Cohen et al. 
2003). On the lower branch, most stars have [Pb/Ba] < 
1.0, except for two stars (CS 22964-161 and HD 198269) for 
which [Pb/Ba] > 1.0. According to Thompson et al. (2008), 
CS 22964-161 is a triple system with a double- lined spectro- 
scopic binary and a third component that might be respon- 
sible for its anomalous abundance pattern. The diverging 
branches of [C/Fc] versus [Fc/H] arc also quite prominent 
in Fig. 8 of Pereira & Drake (2009), who make no com- 
ment on this. Our star G 24-25 belongs to the lower branch. 
Careful inspection of the top panel of Fig. 8 shows evidence 
of the diverging branches in the [/is/Fc] vs. [Fe/H] diagram: 
stars located above the dashed line correspond to the high 
branch in Fig. 9 (except for HD209621, a CEMP-r/s star), 
while those below the dashed line correspond to the low 
branch. Such diverging branches could also occur in the 
[^s/Fe] vs. [Fe/H] diagram, but do not appear to do so, 
as seen in the middle panel, owing to the smaller effect of 
the neutron-exposure and large scatter in the data. This 
similarity seems to be reasonable in the sense that the en- 
hancements of C, Pb, hs, and probably Is all correspond to 
the strength of neutron radiation fiux but to a decreasing 
degree. 

Masseron et al. (2010) completed a holistic study of 
the CEMP stars. They compiled abundances from anal- 
yses of high resolution spectra of 111 CEMP stars and 
21 Ba stars, including both binary and non-binary stars, 
and covering nearly all the chemically peculiar binary stars 
selected in this work. For four stars (HD26, HD 187861, 
HD 196944, and HD 224959), Masseron et al. (2010) give 
updated abundance results. Comparing the left panel of 
Fig. 5 in Masseron et al. (2010) ([Ce/Fe] vs. [Fe/H], where 
Ce is a representative /is-process element) with the top 
panel of Fig. 8 in this work, we can see that [Ce/Fe] shows 
a significant scatter at all metallicities and no obvious de- 
creasing trend with increasing [Fc/H], in contrast to our 
data for [/is/Fe] in Fig. 8. There are two possible rea- 
sons for this difference: the first is a statistical effect where 
[hs/Fe] is based on the average abundance of five s-process 
elements, which thus has a smaller the scatter; the second 
is that Masseron et al. (2010) adopted more stars, i.e. all 
CEMP-s stars (the CEMP-low-s included) in their Fig. 5, 
but we consider only the s-process enhanced binary stars. 

4. Conclusions 

On the basis of high resolution and high S/N spectra, the 
abundances of carbon and 12 neutron-capture elements, in- 
cluding Pb, are obtained for G 24-25. The C and s-process 
element abundances from the first, second, and third peaks 
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show significant enhancements in this peculiar star with re- 
spect to the Sun and also with respect to a reference star 
G 16-20, which has similar atmospheric parameters as G 24- 
25. The abundance ratios [Pb/Fe] = 1.68 and [Pb/Ba] = 
0.33 show that it is not a typical lead star (Van Eck et 
al. 2001), which are often found at low metallicity ([Fe/H] 
< —2) to have [Pb/Ba] > 1.0. Owing to its binary signa- 
ture, high radial velocity, and abundance pattern, we sug- 
gest that G 24-25 is a newly discovered CH subgiant with 
a clear detection of the Pb line at A4057.8. 

Two simple AGB wind-accretion models have been 
adopted to predict the theoretical abundances of the s- 
process elements. The comparison of the observed abun- 
dances with the model predictions indicates that mass 
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transfer via wind accretion across a binary system from 
its principal component, a low-mass AGB star, succeeds 
in reproducing the enhancements of the neutron-capture 
elements. When compared with other chemically peculiar 
binary stars in the literature, we have found that G 24-25 
follows the general trend of [s/Fe], [hs/ls], and [C/Fe] ver- 
sus [Fe/H] established for Ba stars, CH stars, and CEMP-s 
stars implying that there is a similar origin for these ob- 
jects. 
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Table 2. Atomic data and equivalent widths of spectral lines, as well as 
the derived abundances for G 24-25 and G 16-20. 

This table contains the following information: wavelength in angstroms 
(Column 1), element identification (Column 2), excitation potential 
(Column 3), log of the oscillator strength (Column 4), reference for log gf 
(Column 5), measured equivalent width in milli-angstroms and derived 
abundance of G 24-25 (Columns 6 and 7), measured equivalent width in 
milli-angstroms and derived abundance of G 16-20 (Columns 8 and 9). 



Wavelength 


Elem 


E. P. 

(eV) 


log 3/ 


Rcf. 


G 24-25 
EW (mA) 


loge 


G 16-20 




(A) 


EW (mA) 


loge 


4932.049 


C 


7.68 


-1.68 


1 


14.6 


8.08 






5052.167 


C 


7.68 


-1.30 


1 


23.8 


7.97 


1.2: 


6.58: 


5380.337 


C 


7.68 


-1.62 


1 


14.1 


8.02 






4607.340 


Sr 


0.00 


0.28 


2 


18.0 


1.96 


5.7 


1.19 


4077.714 


Sri 


0.00 


0.17 


2 


234.8 


1.79 


154.3 


1.09 


4398.010 


Yi 


0.13 


-1.00 


3 


50.9 


1.58 


18.3 


0.56 


4883.690 


Yi 


1.08 


0.07 


3 


56.1 


1.52 


25.0 


0.60 


4900.110 


Yi 


1.03 


-0.09 


3 


56.2 


1.63 


23.5 


0.67 


5087.430 


Yi 


1.08 


-0.17 


3 


43.4 


1.45 


15.0 


0.54 


5123.220 


Yi 


0.99 


-0.83 


3 


22.0 


1.51 


4.5 


0.53 


5200.420 


Yi 


0.99 


-0.57 


3 


30.8 


1.46 


8.8 


0.58 


5205.730 


Yi 


1.03 


-0.34 


3 


39.3 


1.47 


14.7 


0.64 


5402.780 


Yi 


1.84 


-0.44 


3 


8.6 


1.44 






4687.800 


Zr 


0.73 


0.55 


4 


4.7 


2.14 






4050.330 


Zri 


0.71 


-1.06 


5 


23.2 


2.06 


11.5 


1.47 


4090.510 


Zri 


0.76 


-1.01 


6 


26.7 


2.14 






4161.200 


Zrl 


0.71 


-0.59 


5 


46.7 


2.15 


24.1 


1.39 


4208.990 


Zri 


0.71 


-0.51 


5 


47.0 


2.07 


26.1 


1.36 


4258.050 


Zri 


0.56 


-1.20 


5 


27.2 


2.13 






5350.080 


Zri 


1.83 


-1.24 


6 


2.9 


2.20 






5350.310 


Zri 


1.77 


-1.16 


5 


4.2 


2.22 






5853.688 


Bai 


0.60 


-0.91 


7 


94.6 


2.15 


29.5 


0.55 


6141.727 


Bai 


0.70 


-0.03 


7 


163.1 


2.08 


69.0 


0.53 


6496.908 


Bai 


0.60 


-0.41 


7 


142.5 


2.13 


61.6 


0.64 


4526.110 


Lai 


0.77 


-0.59 


8 


15.0 


1.22 






4662.510 


Lai 


0.00 


-1.24 


8 


22.8 


1.28 






4748.730 


Lai 


0.93 


-0.54 


8 


14.7 


1.23 






4920.960 


Lai 


0.13 


-0.58 


8 


48.8 


1.38 


10.1 


0.06 


4921.780 


Lai 


0.24 


-0.45 


8 


45.2 


1.27 


10.6 


0.07 


4053.490 


Cei 


0.00 


-0.61 


9 


29.5 


1.71 






4062.230 


Cei 


1.37 


0.30 


10 


19.6 


1.88 


1.2: 


0.36: 


4083.220 


Cei 


0.70 


0.27 


9 


36.1 


1.70 






4117.290 


Cei 


0.74 


-0.45 


9 


11.3 


1.69 






4118.140 


Cei 


0.70 


0.13 


9 


32.4 


1.74 


..» 




4120.830 


Cei 


0.32 


-0.37 


9 


31.9 


1.85 






4127.360 


Cei 


0.68 


0.31 


9 


40.9 


1.76 






4153.130 


Cei 


0.23 


-0.80 


9 


16.7 


1.74 






4222.600 


Cei 


0.12 


-0.15 


9 


46.6 


1.80 


6.6 


0.28 


4427.920 


Cei 


0.54 


-0.41 


9 


19.2 


1.71 






4460.230 


Cei 


0.47 


0.28 


10 






10.7 


0.41 


4483.900 


Cei 


0.86 


0.10 


9 


29.0 


1.79 






4523.080 


Cei 


0.52 


-0.08 


9 






3.1 


0.24 


4539.740 


Cei 


0.33 


-0.08 


9 


42.5 


1.78 


5.4 


0.30 


4560.960 


Cei 


0.68 


-0.26 


9 


23.5 


1.82 






4562.360 


Cei 


0.48 


0.21 


9 


50.6 


1.85 


7.0 


0.29 


4565.840 


Cei 


1.09 


0.07 


9 


18.2 


1.74 
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Table 2. continued. 



Wavelength 


Elem. 


E. P. 

(eV) 


log 3/ 


Rcf. 


G 24-25 
EW (mA) 


loge 


G 16-20 




(A) 


EW (mA) 


loge 


4628.160 


Ceil 


0.52 


0.14 


9 


46.1 


1.83 






5044.030 


Ceil 


1.21 


-0.14 


9 


10.4 


1.72 






5187.450 


Ceil 


1.21 


0.17 


9 


18.1 


1.71 






5274.230 


Ceil 


1.04 


0.13 


9 


23.3 


1.72 






5330.540 


Ceil 


0.87 


-0.40 


9 


12.4 


1.72 






5468.380 


Ceil 


1.40 


-0.07 


9 


8.9 


1.75 






5610.246 


Ceil 


1.05 


-0.74 


10 


6.0 


1.85 






5322.800 


Prii 


0.48 


-0.32 


11 


syn 


1.28 


syn < 


-0.38 


4012.700 


Ndii 


0.00 


-0.74 


12 


27.0 


1.56 


4.3 


0.35 


4051.150 


Ndii 


0.38 


-0.30 


13 


28.5 


1.54 


6.4 


0.48 


4232.380 


Ndii 


0.06 


-0.47 


13 


35.1 


1.53 






4542.603 


Ndii 


0.74 


-0.28 


13 


17.2 


1.49 


2.4 


0.34 


4645.770 


Ndii 


0.56 


-0.76 


13 


10.0 


1.49 






5089.830 


Ndii 


0.28 


-1.40 


2 


5.5 


1.53 






5092.780 


Ndii 


0.30 


-0.61 


13 


21.4 


1.46 


2.6 


0.29 


5192.620 


Ndii 


1.14 


0.27 


13 


25.6 


1.53 


3.6 


0.33 


5234.210 


Ndii 


0.55 


-0.51 


13 


21.2 


1.59 






5249.600 


Ndii 


0.98 


0.20 


13 


26.9 


1.47 






5255.510 


Ndii 


0.20 


-0.67 


13 


25.8 


1.53 






5311.480 


Ndii 


0.99 


-0.42 


13 


8.3 


1.44 






5319.820 


Ndii 


0.55 


-0.21 


12 


31.1 


1.55 


3.6 


0.20 


5385.890 


Ndii 


0.74 


-0.82 


12 


7.6 


1.54 






4244.700 


Smii 


0.28 


-0.73 


14 


10.1 


0.93 






4434.320 


Smii 


0.38 


-0.26 


15 


24.6 


1.06 


4.6 


-0.04 


4458.520 


Smii 


0.10 


-0.78 


15 


14.3 


0.97 






4519.630 


Smii 


0.54 


-0.43 


14 


14.8 


1.07 






4523.910 


Smii 


0.43 


-0.58 


14 


13.7 


1.06 






4566.210 


Smii 


0.33 


-0.92 


15 


9.6 


1.11 






4642.230 


Smii 


0.38 


-0.52 


14 


12.1 


0.93 


2.4 


-0.09 


4674.600 


Smii 


0.18 


-0.56 


15 


15.6 


0.90 


3.7 


-0.07 


4687.180 


Smii 


0.04 


-1.17 


15 


9.7 


1.07 






4129.700 


Euii 


0.00 


0.22 


16 


syn 


-0.25 


syn 


-0.50 


6645.130 


Euii 


1.38 


0.20 


16 


syn 


-0.30 : 


syn 


-0.56: 


4191.080 


Gdii 


0.43 


-0.48 


17 


syn 


0.84 


syn 


0.08: 


4057.810 


Pbi 


1.32 


-0.22 


18 


syn 


2.28 


syn 


<0.60 
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(9) Lawler et al. (2009); (10) Palmeri et al. (2000); (11) Mashonkina et al. (2009); 
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